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Abstract. We consider the linear stabilities of the regular n- gon relative equilibria of the 
(1 + n)-body problem. It is shown that there exist at most two kinds of infinitesimal bodies 
arranged alternatively at the vertices of a regular n-gon when n is even, and only one set of 
identical infinitesimal bodies when n is odd. In the case of n even, the regular n-gon relative 
equilibrium is shown to be linearly stable when n ^ 14. In each case of n = 8, 10 and 12, 
linear stability can also be preserved if the ratio of two kinds of masses belongs to an open 
interval. When n is odd, the related conclusion on the linear stability is recalled. 
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1. Introduction 

The famous essay Maxwell(1859)|l| discussed the stability of the rings of Saturn by mathe¬ 
matical modeling. Under the hypothesis that the motion of the rings is uniform and should be 
stable, he abandoned the possibilities that the rings are connected solid, or continuous liquid. 
He concluded that the rings should be composed of countless discrete particles (in. p.66) after 
an important intermediate study. In that study, he conceived a very nearly circular and uniform 
ring of satellites with equal masses, and found that the ring is linearly stable provided the ratio 
of Saturn’s mass to the total mass of n satellites is greater than 0.4352n 2 (m, p.25 or p.59). 
This essay has aroused attentions of many aspects. One can also refer to a French monument 
written by Tisserand[2j. 

Pendse (1935) [3] pointed out that Maxwell’s ring model could not be applied to the cases 
when n is smaller than seven, because Maxwell assumed tacitly the center of primary to be 
the center of masses. Scheeres & Vinh (1991) [4] accepted Pendse’s argument, and analyzed the 
characteristic equations very attentively. They pointed out a miscalculation of Pendse, and con¬ 
cluded that the upper bound of the ratio of the total mass of n satellites to the dominant mass 
for stability can be represented by an asymptotic series of n if n > 7, and is about l/(0.4352n 2 ) 
when n is large. Similar results are also achieved by Roberts (2000) [TD] and Vanderbei & Kole- 
men (2007) [14]. 

This work is supported by the National Natural Science Foundation of China under grant No. 10833001. This 
paper has been published in Qual.Theory Dyn. Syst. (2013) 12: 255-271, DOI 10.1007/sl2346-012-0089-6, and 
is available at http://link.springer.com/article/10.1007/sl2346-012-0089-6 . 



2 


XU 


Consider the planar TV-body problem with N = n + 1 in a uniformly rotating coordinate 
system with the origin at the center of masses. Suppose one mass is large and the other n masses 
are very small. As the order of magnitude of the n small masses tends to zero, the limiting case 
of a relative equilibrium of this problem, is called a relative equilibrium of the (1 + n)-body 
problem, which was defined by G. R. Hall in an unpublished paper [5]. The relationship between 
the relative equilibria of (1 +n)-body problem and those of the ./V-body problem with N = n +1 
has been revealed by Moeckel (1994) [6], and will be mentioned in section 2. 

If all the infinitesimal masses are equal, Casasayas et al. (1994) [7) showed that when 
n > e 73 , the only stationary configuration is Maxwell’s ring configuration. When n < 8, other 
configurations exist, referring to Salo & Yoder (1988)[8] . According to Salo & Yoder’s numerical 
quest, it seems that there is only one type of stationary configuration when n > 9, and this 
conjecture is also supported by the numerical experiment of Carles Simo (0, p.326). When 
2 < n < 4, the numerical results on the numbers of stationary configurations also concide 
with analytical proofs, see Cors et al. (2004) [D], Albouy & Fu (2009) [121 ■ For the cases of 
n = 5,6, 7,8, 9, they seem to be much more difficult to deal with, as the case of n = 4 had 
already been claimed by Hall to be difficult to handle ( 0 , P-12). 

If the infinitesimal masses are arbitrarily given, Maxwell said that “we must calculate 
the disturbing forces due to any given displacement of the ring” (ffl, P-38), so it is difficult to 
determine the stability of a non-regular configuration. Salo & Yoder (1988) studied a special case 
of n = 3 with a background of Saturn’s coorbital satellites, and Renner & Sicardy (2004) pT] 
studied the cases of n = 3,4,5. Recently, Corbera et al.(2011) [13] researched into the case of 
n = 3 and obtained two classes of new configurations generated by changing the infinitesimal 
masses. They also calculated the number of configurations with a high precision, and found 
that the number varies from five to seven. 

The new results gotten in this paper are especially about the existence and linear stability 
of the regular n -gon relative equilibria in the (l+n)-body problem with two kinds of infinitesimal 
masses. In section 2.1, we introduce Hall’s potential function on the regular n-gon relative 
equilibria of the (1 + n)-body problem. Then the method for calculating eigenvalues of the 
Hessian matrix of Hall’s potential function is described in section 2.2 and 2.3. In section 3, 
we show the existence of a set of positive infinitesimal masses for the regular n-gon relative 
equilibria. In section 4, we show the linear stabilities of regular n-gon relative equilibria with n 
even. The last section is the conclusion part, and we also give some vistas for future study. 


2. Regular n-gon relative equilibria 

There are two ways for the definition of a relative equilibrium of the (l-t-n)-body problem. One 
way is to make a first approximation of the full system, then write out the conditional equations 
for a relative equilibrium, for example, Salo, Yoder (1988). The other way is to make a limiting 
case of the conditional equations for a relative equilibrium of the full system, for example, 
Moeckel (1994). We obtain the conditional equations for relative equilibria of the (1 +n)-body 
problem by the second way. 

2.1. Relative equilibria 

Consider the Newtonian iV-body problem with N = n + 1 in a rotating coordinate frame 
with the origin at the center of masses. Set the angular velocity of this frame as a constant 
ui = 1. Suppose there is only one dominant mass and the other n masses are very small. Denote 
the dominant mass as mo, and the other masses as TOj, i = 1, • • • , n. Introduce the positions 
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and conjugate momentums of the small masses as Xi,yi £ R 2 , i = 1, • • • 
Hamiltonian function can be written, 


where 


n 


i= 0 




n— 1 n 


EE 

2=0 j=i-\-l 


miirij 

I \xj -ar*ir 


n 

m 0 = 1, mi = £fj,i , #o = ~y 

2=1 


n 

2/o = “E^ 1 

i=l 



and the prime represents transposition. 


, n, separately. The 


( 2 . 1 ) 


0 

-1 


1 

0 


Let’s recall a theorem given by Moeckel (1994) in order to understand the relationship 
between the relative equilibria of the (1 + n)-body problem and those of the iV-body problem 
with N = n + 1. 


Theorem 2.1 (Moeckel, 1994). Let x e be a family of relative equilibria of the N-body problem 
with N = n + 1, too = 1, and irii = £fii, 1 < i < n. Suppose that as e —> 0, x e converges to a 
non-degenerate relative equilibrium, x, of the (1 + n)-body problem. Then x e is non-degenerate 
for s sufficiently small. In this case, x e is linearly stable for e sufficiently small if and only if x 
is a local minimum o/V, which is defined in (| 2. below. 


The following proposition can be found in both Hall (preprint) and Moeckel (1994). 
Proposition 2.2. For x as given in the theorem above, we have xq = 0, \xt\ = 1, i = 1,2, • • • ,n. 

The conditional equations for a relative equilibrium of Hamiltonian system can be 
written, 


. _ &H _ Vi r _ n 
Xi — „ — + J Xi — 0, 

oy% m% 

. _ _ J mjmo(x o - Xi) 

dxi Vl ||x 0 — Xi \\ 3 -“t- b. || Xk-Xi 

" " «=1,^2 


mim k {xk - xQ _ Q 


( 2 . 2 ) 

(2.3) 


where i = 1,2, ■ ■ • , n. 

One gets Jyt = iriiXi from (12.21) . and substitutes it into (12.31) . Let Xi = ri (cosi?*,sin!?,;) / , 
(12.31) make inner product with (— sin'!?,, cosi?i) / , and e — > 0. Then the conditional equations for 
a relative equilibrium of (1 + n)-body problem can be obtained, 


E Tk sin(i?fc - !?i) ( ^3-1 J = o, 

k=l,k^i 


(2.4) 


k.i 


where r\-,i = 2 sin ^ 2 ^ ^ , k ^ i. 

Integrate the left hand side of equation (ES of di, yields 

1 rl.\ 


n— 1 n 

v = E 

i=l fe=i+l 


Tk,i 


(2.5) 


and function V is called Hall’s potential function. 
Equations (12.41) refers to a function below, 

1 


IFfp) = sin I 1 — 


8| skr tpj 2| 


ip £ {</?|<^ 2kn, k £ Z}. 


( 2 . 6 ) 
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The first order derivative of fF(<p) is 

f (y) =?'&) = oi 1 y| ( ~- 1 ) + cos ip. (2.7) 

81 sin f\ y sm f J 

Note that functions f((p) and fF{ip) are useful for later analysis, so let’s recall a proposition, 
which can be found in 0,0,0,1123- 

Proposition 2.3. Functions fF(ip) and {(ip), with the field of definition {ip\<p £ K, 2kir,k £ 
Z}, satisfy 

1 . .F(-v>) = -^(^), ^ ± 2 n) = ?(<p), F(n - ip) = -F(ip - n), 

2 - f(~¥>) = %), f(7T-<^) = f(v3-7r), %)>f(7r) = -|, 

3. 3 ci > 0, such that {(<p) > Ci/ip 3 — 1, for ip £ (0, 7 r), 

4. f"M=>"(^)>0, V^£ (0, 27 r) 

These properties of the two functions are shown in the figure CD 

2.2. Circulant matrix 

The linear stabilities of the relative equilibria of the (1 + n)-body problem depend on the 
eigenvalues of the Hessian matrix V##. When all the infinitesimal masses are equal, the matrix 
is a circulant matrix, and Moeckel (1994) concluded that such regular n-gon relative equilibria 
are linearly stable if and only if n > 7. The definition of the circulant matrix is given below. 

Definition 2.4. Let j be a positive integer. Matrix A of j x j is a circulant matrix, if it satisfies 

Ai,k = ^4{j +i}j,{k+i}ji i,k,l £ N, 

where is an element of the matrix and {i + l}j represents i + l modulo j. 

In a regular n-gon relative equilibrium, suppose one infinitesimal body lies on the positive 
x axis, and count them from this one by the anti clockwise direction. The angle positions of the 
infinitesimal bodies are 

fl« = 1)77 , t = 1,2,(2.8) 

n 

Let’s review the following proposition 12.51 about the circulant matrix, and this propostion 

can be found in [123- 


Figure 1 


F(<p) represents function (12.61) and f (ip) represents function (12.71) 
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Proposition 2.5. Every circulant matrix A with rank j has eigenvectors 

£(0 = -L (l, e -2^i(i- 1 )/3 } ... , e -2«-i(i-i)0--i)/j-y } i = i ; 2} ... , j, 

where i = >/— 1. 

Denote the element in the i-th row and k-th column of matrix A as Aj ) / c . The eigenvalues 


ai=53 A 1 .fc e_2,ri(, " 1)(fc " 1)/j ’> ^ = !,2,- • • , j. 

fc=i 

One has A = Udiag{cq}U _1 , where U has the eigenvectors as columns in order, and diagjcq} 
is a diagonal matrix of eigenvalues of matrix A. 

Let B 6e another circulant j x j matrix, and denote its element in the i-th row and k-th 
column as Bif.. Its eigenvalues are 

Pi=i2 B Z = 1,2,---,i, 

k=1 

respectively. Then 

1. A and B commute and AB = BA = Udiag{az/3j}U _1 , AB is also a circulant matrix. 

2. A + B is a circulant matrix and A + B = Udiagjcq + /3;}U 1 . 

3. If ai 0, l = 1,2, • • • , j, then A is nonsingular and 

A- 1 =Udiag{cq- 1 }ir 1 . 

However, if not all the infinitesimal masses are identical, will no longer be a circulant 
matrix. 


2.3. Eigenvalues of a block circulant matrix 

Define a block circulant matrix S as below, 

S = 


A C 
C T B 


(2.9) 


|AI 2i - S| = 


where A, B and C are all circulant matrices of j x j. The upper T represents transposition. 

Denote I- 2 j as a 2 j x 2 j unit matrix, and I as a j x j unit matrix. The characteristic 
polynomial of matrix S can be calculated by applying some elementary transformations. 

AI — A C 

-C T AI — B 

I 0 

— C r (AI — A) -1 I 

AI-A 0 

0 (AI-B)-C T (AI-A)-!C 

= |(AI — A) [(AI-B)-C t (AI-A)- 1 C]| 

= |A 2 I - (A + B)A + [AB - (AI - A)C T (AI - A) -1 C] |. (2.10) 

The eigenvalues of matrix C are 


AI - A -C 

0 (AI-B)-C T (AI-A)-!C 

I —(AI — A) -1 C 

0 I 


7 i = £c 1>fc e- 2 ’ ri(, - 1)(fc - 1)/ J 1 1 = 1,2 ,-•• ,j, 


k—1 
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where C\ t k is the element in the first row and the k -th column of matrix C. The following 
equation can be obtained by Proposition 12.51 

U -1 (AI - A)C t (AI - A) _1 CU = diag{7;7;}, 

where 7 represents the eigenvalues of matrix C T , and equals to the conjugate of 7. So, charac¬ 
teristic equations can be obtained as below, 

Xf - (ai + Pi)\i + aifii - 7;7z = 0, 1 <l<j. (2-11) 

One can also find that 


ai — atj + 2~i, Pi — Pj+2-u 111 1 — lj+2-ilj+2-i- 

If cq, Pi and iih are calculated for 1 < l < j, then the eigenvalues of the block circulant 
matrix S can be calculated by solving equations (12.111) . 


3. The set of masses for regular rt-gon relative equilibria 


Consider a regular n-gon relative equilibrium. There is only one set of identical infinitesimal 
masses for a regular n-gon relative equilibrium when n is odd, and there exist two kinds of 
infinitesimal bodies arranged alternatively at the vertices of the regular n-gon when n is even. 
These can be shown to be true by referring to a proposition given by Renner & Sicardy (p.403) 
and the method in section 2.3. 


Let’s review that proposition in this paragraph. With the help of equations (12.41) and (12.61) . 
one can write conditional equations for the relative equilibria of the (1 + n)-body problem in 
the following way, 

n 

y hkF{dk - Pi) = 0 , i = 1,2, (3.1) 

k=l,^i 


Suppose that angles ($ 1 , • • •, i3 n ) correspond to a relative equilibrium. Define ~ $i) 


with k 

7^ b 

kk 

= 0 , 

and note that 

matrix 

form, 






( 

0 

.F 12 

^13 




T 12 

0 

F 23 



\ - 

In 

2n 

•F3 n 


2 n 

0 / 


\ 

( hi ^ 


( hi \ 


h2 


h2 



= M n 


J 

V hn ) 


V hn J 


= 0 


nx 1 • 


(3.2) 


If n is odd, and rank(M ra ) is n — k, where n is an odd integer, then there exists a K-parameter 
family of mass vectors (/xi, • • •, If n is even, then rank(M„) is generally n, and there is 
generally no family of mass vectors for a relative equilibrium. However, Renner & Sicardy just 
mentioned examples of when n = 2,3 in their remarks. 


Substitute the angle positions given in (12.81) into (13.21) . For a regular n-gon relative equi¬ 
librium with n > 3, we will show that, 


Proposition 3.1. When n is odd, rank(M„)= n—1, then there exists only one parameter for the 
set of infinitesimal masses. When n is even, rank(M n )= n—2, then there exists a two-parameter 
family of mass vectors. In addition, the two kinds of masses are alternatively arranged at the 
vertices. 
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Proof. The rank of matrix M„ equals to the number of non-zero eigenvalues. For a regular 
n-gon relative equilibrium solution, one finds that 

fa = - 00 = -H$i - #i) = -fji, 

Pik Pi+l mod n, k+l mod n • 

The two expressions above verify that the matrix M n is both asymmetric and circulant. Note 
that 


27T 


Pi k + Pi, n + 2 -k = P — (k - 1) + P[ —(n + 1 - k) = 0. 


2n 


The eigenvalues are 


^ ' P\ k * ^ 


. P ~ 2iri(J-l)(fe-l)/n 


fe=1 

1 


= -^2 \F lk ■ e -27ri(i-l)(fe-l)/n + jr l>n+2 _ fc . e -2ri(M)( n +l-fc)/ n 
Z fc =2 
1 " 

= 2^^' 

k—2 

. \r-' -j-, 2A’7t . 2(/ — l)/c7T 


g—27ri(Z— l)(fc—l)/n 27ri(Z —1)(1 — k)/n 


k =1 


n— 1 


- E 


I sin ■ 

n ' n 


2lkir (4 — 2l)k7r 
cos-cos- 


fc=i 


1 ^-4 cos ^ . 2(1 - l)kn 


-Y 


4 f- sin 2 ^ 

fc=1 n 


sin ■ 


= — i/i(n, Z), 1 < l < n. 


One finds Ai = 0 and A n m+i = 0 when n is even. It is evident that A; = —A„+ 2 -; when 
2 < l < n, so one can just consider 2 < l < [(n + 1)/2J. Here |_-J represents the biggest integer 
that is no greater than the number inside the symbol. 

The first part of f\(n,l) can be simplified as 


1 

2 


n— 1 


E 


2lkn (4 — 2l)kn 
cos-cos- 


if/= 2, 

0, if / = 3,4, ■ • • , [(n + 1)/2J, 


and the second part of fi(n, l) can be denoted as / 2 (n, / ), 


h(n,l) 


n— 1 

E 

k= 1 


COS — 


sin 


2 /C7T 


sin 


2(1 — l)fc7T 


n 


n> 3. 


To check the convex property of the function f 2 (n, /), one may judge whether the following 
expression is positive or negative. 

_ L"jr /I /1 1 7T 

f 2 (n, l) + f 2 (n, 1 + 2)- 2f 2 (n, l + 1) = - V cos — sin —, 

z ' n n 
fc=l 

where 2 < l < [(n — 3)/2j, n> 7 . 

One can find that the above equality is smaller than zero, because, in Yk=l cos ~ si n 
the positive part is greater than the absolute value of the negative part. So function f 2 (n, l) is 
convex with the mouth downward. 
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The two ends of /i(n, l) can be shown to be positive when n is large. When n is odd, 
is an integer, and 


n + 1 l ^cos^f . n-l i 1 (-l ) fe+1 

£ —^ >0 - 


4 ^ sin 2 — n 

k= 1 n 


2 ^ tani u 


When n is even, = §, and 


n, 1 cos . n — 2 




■ sin - 


2' 4 ^ sin 2 ^ n 

fc=l n 


n/2_1 cos 2 ^ 

**= £ —^(-i) fe+ 1 > 0 . 


fc= 1 


sin — 


Thus fi(n, L(n + 1)/2J) > 0 when n > 3. 

On the other hand, we’ll check whether /i(n, 2) is greater than zero when n is large. 


fi(n,2) 


n— 1 

n 1 v-^ 

-■>+?£ 


fc= 1 



Apparently, 1/ sin a’ is convex with the mouth upward as a 6 (0, 7 r), so one can have the following 
relationship by the trapezoidal rule, 


n— 1 1 

x— 

J cm . 


fc=i 


7r 

. - > 
sin — n 


> 


ii ^ i 


2 sin 

(n-l) 


fc =2 


sin 


1 

hi. + 2 


smx 


-dr = 2 In ■ 


sin ■ 

cos ii 
sm ^ 


1 — 1)7 r 


7T 

n 


One can also get the summation of the series below, 


Tl — 1 j -i , jr 

. &7T 1 + COS 

sm — =-;——— 

n sm ■ 


k= 1 


n 


sm 


2 n 


Define a function as below, 


7 r 1 cos a: cos a; 

= + ln “ 2 sin a; 

7 T 1 , COS X 1 

>- 1 -In- 

4a: 2a’ sin x 2 sin x 


x 


1 r, cos x 

— In- 

2x L sin a sin a 


(3.3) 


(3.4) 


The definition domain of / 3 (a) is (0, 7 t/2 ), and the zero root depends on the function / 4 (a) as 
below, 


/ 4 (a) = In 


cos a 
sin a 


2 a 7 r 

sin a 2 


Let a = 7 t /2 n, and function fiin) is a monotonously increasing function as n > 3. One can find 
/4(41) < 0 and /4(42) > 0, then one has /i(n, 2) > 0 if n > 42. 


In addition, one can use one calculating software to calculate and find that there are no 
very near zero values for fi(n,l) when 3 < n < 41 and 2 < l < |_(n + 1)/2J. So, there is only 
one zero eigenvalue Ai if n is odd, and there are only two zero eigenvalues when n is even. This 
completes the proof. □ 
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4. Linear stabilities of the regular n-gon relative equilibria when n is even 

Exact expressions of the elements of the Hessian matrix V## are given as below, 

3 + cos($z - i?j) 


V'di'dl — fj'ifj'l 


2 rL 


+ COs(l?z - 1?i) 


l ± i, 


2 j 

^ Vdi$ v 
1 = 1 , 


Especially, note that = —HiHi^i^i — $i). Consider that there are just two kinds of infini¬ 

tesimal masses, let Hk = Hi if k is odd, Hk = H? if & i s even. 

V #,9 is equivalent to the block circulant matrix S in (12.91) by row and column exchanges. 
The elements in the four blocks of S are 

-Ai^k — = l / i?2i^2fc) Ci,fc = f / l?2»-H?2 k> 

and 1 < i, k < j. Matrix S is a symmetric matrix, so all its eigenvalues are real. Also note 
that ai , pi and qzy; are real. The aime of this section is to show that there exist no negative 
eigenvalues in S such that the linear staibilites of the regular n-gon relative equilibria can be 
assured by Theorem 12. II 

The sufficient conditions for all the roots of equations (12.111) to be nonnegative are 

otl+ Pl> 0, aipi > yzy i = |7| 2 , 1 <l<j. (4.1) 

These conditions above will be checked by applying the Proposition 12.51 
The expressions for cq, Pi, 7 1 and 7 1 in (14.11) are given, 

ai = Vw, + J2 Vtf ltfM _ 1 e- 2 ’ ri(I - 1)(fc - 1)/ ^, 
k =2 

Pi = V» 2 0 2 + V^e- 2 ^ 1 - 1 ^- 1 ^, 

k =2 


7 1 ='52v.0 1 0 2k e 2m(i 1)(fe 1)/j , 7i = y^2k- I e 


-27ri(Z-l)(fc-l)/j 


k= 1 


k= 1 


where 


v#,#, = hi f f 

fc =1 \ 3 J fc=i ' 3 


2 fc- 1 


2 v—' f 2kn\ \r^ r 

V# 2 tf 2 = 1*2 2 ^ f l ~ f 


fc=i 


fe=i 


2 fc — 1 

j 


Note that a; and /?; can be separated into two parts as described below, 

on = Hi SiCb 0 + HiH2g2(j), 

Pi = HlgiUJ) +HiH2g2{j), 
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where 


k =1 ' J 

3 


3 


g 2 (i) = 


fc=l 


2 fc — 1 


-7T 


J 


The exact forms of 7 ; and 7 ; are also given, 


7i = -/ii/i 2 X! f 
k=l 
j 

7 1 = -M1M2 X] f 

fc=l 


2k - 1 
- 7 

j 


2k - 1 


3 


2TTj(l — l)(fc — 1) 




7 T e 


The expression of | 7 ;| can be calculated in the following way, 


7*7* = ( 2 4^) 2 - ( 2 V 1)2) 


and one gets 
where 


1721 = MlM2g3(j,0> 


gsO ',0 = 


fe=i 


2 /c — 1 


J 


7T cos ■ 


(Z — l)( 2 fc — 1 ) 


The following expression of x(j, l) is also given as it is necessary and crucial, 

X(j, 0 = - 7;7 j ) 

M 1 M 2 

= M 1 M 2 [g+ g l(j) - g§(i, 0 ] + (m? + M 2 )gi( j) 0 g 2 (j)• 


(4.2) 


Some propositions about gj (j, /), g 2 (j ) and g 3 (j, Z) are given such that one can check the 
inequalities in (Ell- 

Proposition 4.1. gi(j, l) increases with l when 2 < l < \_(j + 1 )/ 2 J and j > 3, where |_-J means 
the maximum integer which is no greater than the real number inside. 


Proof. 


where 


3 -1 


2/C7T 


giO'iO = 


fc =1 


3 + cos —j- 2kn 

- ~-r -b cos- 

16 sin 3 — j 


1 — cos 


2(1 — l)fc7T 


= ci(i) + c 2 (0, 


i 3 + cos 


2kn 

j 


Cl ® ^ 16 sin 3 ^ 

fc=l J 


1 — cos 


2(1 — l)kn 


3 -1 


2&7T 


C2 


(o = E c °s— 1 - 


fe=i 

j 


cos 


E 2 k / 

cos 1 _ 


fc=l 


cos 


2(1 — 1)/C7T 


2(1 — 1)/C7T 


— 


COS 


2(1 - 2)kn 


k= 1 
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It is easy to find that 02 ( 2 ) = — and C 2 (l) = 0 if l 2. Apparently ci(Z) is positive. In 
addition, 


Cl 


(Z + 1)- C1 (0 = ]T 


J 1 2 - sin 2 !f ^ (21 - ljfcir 


, , 4 sin 2 ^ 
fc=i r 


■ sin ■ 


J 


> 0 , 


so ci(Z) increases with l when 2 < l < \_(j + 1)/2J. The conclusion is thus apparent. 
Proposition 4.2. For all j > 7, we have gi(j, 2) > 0. 


□ 


Proof. 




J 


One refers to the relationships of (13.3D and (13.41) . and finds 


7 7 COS 7T- 1 cos :A 

gi(i,2)>-U^-ln-4-i^4 

2 7T sin 4 sin 


Let x = and one has 


2 a 


. _ . cos x x cos x 

hi(x) = — I In —-—- 

2x V sm a; 2smi 


2 a 


“ (l) = h ( 

— ( 

2x V 


cosx 


X 


h 2 (x) = In 


sin x 2 sin x 
cos x X 7T 


■f) 

f). 


sin x 2 sin x 


h^{x) is a monotonously decreasing function of x € (0, n/2), so gi(j, 2) is an increasing function 
of j > 2. One finds gi(6, 2) < 0, and gi(7, 2) > 0. So gi(j, 2) > 0 for all j >7. □ 

Proposition 4.3. For all j > 1, g 2 (j) > 0. 


Proof. As 


S2 (j) 


J 

' 4 sin ^ 
t=l 4Sm 2j 




7 r sin 2 2 

2 a 


--)>0, if j > 1. 


□ 


Proposition 4.4. For all j > 2, g|(j) - g§(j, 2) > 0. 


Proof. One has 


gil?) - g§0') 2 ) = [g2(i) + g3 {j, 2)] • [g2(j) - g3 (j, 2)], 


g2(j) + gs (j, 2) = f (~7“ 7r 




fe=i 


sin 2 2 ^ 17 T 


2k- 1 

1 + COS-;-71 


— 1 


sin 2 ^. 1 7 r N 


2 a 


2 sin ^L_!- 7 r 

. 2 a 


> 0 , 
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and 


... J^ c f2k-l \ ( 2k-1 \ 

g 2 (j)-g 3 (j, 2 ) = (-t—ttJ f 1 cos— - —ttJ 

fc=l \ 3 ' \ J / 

= v ( 1 _ sini ^r 7T _ A 

£A 2sin ^r" 4 V 

>im^i_I_ i. 

7 r sin -fr 4 sin ff, 2 

The function 

1 / cos | x 7 r\ 

monotonously decreases as x £ (0, 7 r). So g 2 (j) — g 3 (j, 2) increases with j if x = jj. As g2(2) — 
g 3 (2,2) = 3^1 - 1 > 0. So, g 2 (j) ± g 3 (j, 2) > 0 when j > 2. □ 


As a result of the above propositions, xti, 2) > 0 is shown when j > 7. Note that when 
l = 1, gi (j, 1) = 0, but one can still find that x(j, 1) > 0 because of Proposition 14.41 Also note 
that when 2 < l < [(j + 1)/2J and j > 3, gi(j, l ) increases with Z, so x(j, 0 increases with l as 
well. In all, one has x(ji 0 > 0 when j > 7. 

If 2 < j < 6 and one still looks for linear stability conditions, one just needs to check 
whether x(j, 2) > 0 and a 2 + /?2 > 0. In these cases, one gets gi(j, 2) < 0, and the discriminant 
of Xti, 2 ) is 

hi= [g?(i, 2) -h g|(j) — g§(j, 2 )] 2 — 4 [gi(j, 2)g 2 (j )] 2 - (4.3) 

/14 is calculated to be positive when 4 < j < 6 , and negative when 2 < j < 3. So when 2 < j < 3, 
Xti, 2) < 0. 

When 4 < j < 6 , the intervals of ^i//x 2 can be obtained by calculating the roots of the 
equation xti, 2) = 0 easily. Such intervals can ensure the correctness of the inequalities in (14.111 . 
These conclusions are summarized as below, 


Theorem 4.5. Given a regular n-gon relative equilibrium of the (1 + n)-body problem with n = 2j 
and j > 2. If there are two kinds of infinitesimal masses pi and /x 2 , the two kinds of masses 
should be distributed alternatively. It is shown that such a configuration is linearly stable if and 
only if j > 7, or one of the following conditions should be satisfied, 


AT / —Its + v^4 —h5 — \fhl \ ■ — 4 5 f, 

M 2 \2gi(j,2)g 2 (j)’ 2gi(j,2)g 2 (j)y 3 

where is given in and 

hs = slti, 2) + g 2 (j) - slti, 2 )- 

The intervals in o can be calculated as 

{ (0.39601454048825,2.525159805412902), if j = 4, 

(0.16709497914366,5.984620274797297), if j = 5, 

(0.061964963348688,16.13815204525851), if j = 6 . 


(4.4) 














Linear stability of the n -gon relative equilibria of the (1 + n)-body problem 
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5. Conclusion and Discussion 

We’ve considered the linear stabilities of the regular n-gon relative equilibria of the (l + n)-body 
problem. When n is odd, we’ve shown that there is only one kind of infinitesimal mass. We recall 
the conclusion of Moeckel(1994) that the regular n-gon relative equilibria with all infinitesimal 
masses equal are linearly stable if and only if n > 7. We have also shown that, if n is even, there 
can exist two kinds of infinitesimal bodies arranged alternatively at the vertices of the regular 
n-gon. Such relative equilibria with n even are linearly stable if and only if n > 14, or the ratios 
of the two kinds of masses must be in given intervals for n = 8,10,12. See the Proposition 13. II 
and the theorem 14.51 

The allure of the relative equilibria of the (1 + n)-body problem comes from both mathe¬ 
matical and astronomical interests. The features of stationary configurations of the (l + n)-body 
problem are fascinating, but it is difficult to explain those features. One question is what they 
are like and how many there are for n given and not necessarily all equal infinitesimal masses. 
Only special cases can be analytically solved for this question, as the stationary solutions are 
generally difficult to calculate. The second one is the inverse one, which is to look for a set of 
masses for n given positions on a circle. 

There are also some questions about the stability. The effects of eccentricities of the orbits 
of the infinitesimal bodies and the oblateness of the primary can also be considered on the 
linear stability of the relative equilibria of the (1 + n)-body problem. 

In addition, as a result of the relationship between the relative equilibria of the (l+n)-body 
problem and those of the TV-body problem with N = n + 1, we can consider the existence and 
stability of some periodic solutions near the relative equilibria of the iV-body problem with two 
kinds of small masses, by referring to Meyer & Schmidt (1993)[15; and Pascual (1998)[16]. The 
nonlinear stability of the relative equilibria or periodic perturbations of the iV-body problem 
with N = n T 1 is also deserved to be researched, and one can try to analyze the properties of 
higher orders or do some simulations by using a symplectic integrator. 
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